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In this paper, we propose the quantum spin pumping in quantum spin systems where an ap- 
phed electric field {E) and magnetic field {H) cause a finite spin gap to its critical ground state. 
When these systems are subject to alternating electromangetic fields; {E, H) — (sin cos 
and travel along the loop Fioop which encloses their critical ground state in this E-H phase 
diagram, the locking potential in the sine-Gordon model slides and changes its minimum. As a 
result, the phase operator acquires 2tt holonomy during one cycle along Fioop, which means that 
the quantized spin current has been transported through the bulk systems during this adiabatic 
process. The relevance to real systems such as Cu-benzoate and Yb4As3 is also discussed. 

KEYWORDS: field induced spin gap systems, quantized spin transport, Dirac monopole, sine Gordon 
tlieory, phase operator, topological stability, Landau Zener tunneling 
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1. Introduction 

Transport phenomena in magnets such as the colos- 
sal magnetoresistance^ and anomalous Hall effect^"^ 
are long-standing subjects in sohd state physics. Re- 
cently, exotic characters of the anomalous Hall effect 
are closed up in ferromagnets,^"^^ antiferromagnets^'^ 
and spin glass^^' with non-coplanar spin configura- 
tions and coUinear ferromagnets with strong spin-orbit 
couplings. ^^^^^ Their common origin is gradually rec- 
ognized as the topological character of magnetic Bloch 
wavefunctions in its ordered phase, sharing the same 
physics with the two dimensional (2D) quantum Hall 
effect. Thereby the spontaneous Hall conductivity is re- 
lated to the fictitious magnetic field defined in the crystal 
momentum space and its quantized value is directly re- 
lated to the first Chern number associated with the fiber 
bundle whose base space is spanned by crystal momenta. 

Meanwhile, in early 80's, Thouless and Niu proposed 
the quantized adiabatic particle transport (QAPT) in 
gapped fermi systems, where an adiabatic sliding motion 
of the periodic electrostatic potential in one dimensional 
(ID) system pumps up an integer number of electrons 
per cycle. ^^'^^ In fact, this quantized particle transport 
is another physical manifestation of the topological char- 
acter associated with the Bloch wavefunctions. Specifi- 
cally this particle (polarization) current is related to the 
fictitious magnetic field defined in the generalized crys- 
tal momentum space which is in turn spanned by the 
crystal momentum and deformed parameters. Because of 
this analogy to the 2D quantum Hall currents, this ID 
quantized particle current is known to be topologically 
protected against other perturbations such as disorders 
and electron-electron correlations.^^ Due to this pecu- 
liar topological protection, the QAPT became recently 
highlighted in the realm of the spintronics, where elec- 
tron pumpings have been experimentally realized in meso 
and/or nano-scale systems. ^^"^^ However the quantized 
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value observed in these experiments should be attributed 
to the Coulomb blockade^^'^^'^^ and is different form the 
aforementioned Thouless quantization, which is assured 
only in the thermodynamic limit. ^'^ 

In this paper, based on the original idea of the QAPT, 
we propose the quantized spin transport (quantum spin 
pump) in the macro-scale magnets. There we discuss sys- 
tematically the topological stability of its quantization. 
Specifically, we clarify in this paper the relation between 
its quantization and holonomy of the phase operator of 
sine-Gordon theory in ID quantum field theory. By us- 
ing this relation, its topological stability can be quanti- 
tatively judged by seeing whether the expectation value 
of this phase operator acquires 27r holonomy during the 
cyclic process or not. 

As for the experimental relevant systems to our theory 
of quantized spin transports, we have the S = ^ quan- 
tum spin chain such as Cu-benzoate and Yb4As3 (charge 
ordered phase) , where its unit cell contains two crystallo- 
graphically inequivalent sites. Accurately speaking, both 



the translational symmetry by one site (T : Sj 



'i+i 



and the bond-centered inversion symmetry which ex- 
changes the nearest neighboring sites (/bond : ^i-j ^ 
Si+j-i-i are crystallographically broken. As a result of 
this peculiar crystal symmetry, the g-tensors of its two 
sublattices are in general different in these systems; 

N 



H — J ^ Si ■ Si+i + 



H-[[g^] + {-mg,]]-S., (J>0)(I) 



where [gu\ and [ga\ represent the uniform and staggered 
component of the g-tensors respectively.^^ Late 90's, 
these quantum spin systems, whose ground states are 
nearly critical because of their strong one dimensionali- 
ties, are experimentally revealed to show the spin gap be- 
haviors under the external magnetic field. Furthermore, 
following theoretical works showed that the field-induced 
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gap behaviors (A ~ H'^/^) are originated from this stag- 
gered component of the £f-tensor, where the effective stag- 
gered magnetic field in the eq. (1) endows its critical 
ground state with a finite spin gap.^^'^" 

On the other hand, the exchange interaction J in 
these S = 1/2 quantum spin chains, JJ2i'^i ' '^'i+ii 
does not have an alternating component, since these 
systems are invariant under the site-centered inversion 
symmetry which exchanges the nearest neighbor bonds 
(Jsiic : Si^j ^ Si+j). However, when we break this sym- 
metry by applying an electric field E along an appropri- 
ate direction,^^ the exchange interaction in general does 
acquire a staggered component; 

N 
i=l 

Futhermore a site-centered inversion operation with the 
sign of E reversed requires that A must be an odd func- 
tion of E. This dimerizing field A is also expected to 
induce a finite mass to its critical ground state, causing 
the spin-Peierls state. 

Based on these observations, we propose a method of 
generating the quantized spin current in this type of spin 
systems (T, /bond : broken, igite : unbroken) by using 
electromagnetic fields. Specifically, we study a following 
5 = 1/2 Heisenberg model with time-dependent stag- 
gered field and bond alternation: 

N 

H{t) = ^[j-(-irA(i)]s,.Si+i 

i=l 

N 

+ h,,{t)j2{-iys!. (3) 

i=l 

where staggered Zeeman field hst{t) and bond alterna- 
tions A{t) are supposed to be controlled by the applied 
electromagnetic fields. In § 2, we will argue by using 
bosonization technique that the quantized number of z- 
component of spins are transported from one end of the 
system to the other during one cycle along the loop which 
encloses the critical ground state at (A, hgt) = (0, 0) in 
the A-hst plane. 

In order to uphold this bosonization argument, we also 
demonstrate the mimcrical calculations in § 3, where we 
evolve the ground state wavefunction according as the 
above time-dependent Hamiltonian. Thereby we attibute 
the quantized spin transport to the Landau-Zener tun- 
nelings^° wliicli indeed happen between several energy 
levels during this cycle process. 

In addition to eq.(3), we also have a uniform compo- 
nent of Zeeman fields in real systems, 

N 

H'{t) = J2H{t)-Sj, 

i=i 

whose magnitude is usually larger than that of the stag- 
gered component /ist(i)- This uniform field, especially its 
X 01 y component, seems to easily flip the spins accu- 
mulated at both boundaries and might spoil the physical 
consequence of the spin transport. Then we also discuss 
in § 3 the effect of this uniform field and the remedy 
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against it, where it turns out that the sweeping velocity 
should be appropriately chosen so as to avoid the relax- 
ation process via this uniform Zeeman field. 

2. Bosonization 

For clarity, let us first illustrate the physics of our 
quantum spin pumping in a simple limiting case; 

= -i:Y.^StSr-_^, + C.C.), (4) 
^ i=i 

//dim = -^Y.^-ir{stS-^^ + S-St+^), (5) 

i=l 

i/st = K,Y.^-ir^t, (6) 

i 

where we take the direction of the staggered Zeeman field 
as the z-direction.^^ In this simplification, we neglected 
the exchange coupling between the 2:-component of spins. 
In this limit, we can get a quadratic form of Hamiltonian 
in terms of the Jordan- Wigner (JW) fermion Sf = flfi — 

\,S+ = //exp(i7rE-=i/l/,)-''''^ Then Fxy forms a 
cosine band in its momentum space, e{k) = — Jcosfca 
{k,a are crystal momentum and lattice constant). In the 
absence of A and /i^t, the fermi points locate at fc = ±^ 
since fermions in the ground state fill up all the k points 
but those with positive energy e{k). 

Nonzero A and/or hgt introduce a finite gap at these 
two Fermi points. The dimer state induced by a finite A 
can be understood as a Peierls insulator^ where the JW 
fermions occupy the bonding orbitals between the two 
neighboring sites and form a valence band in its fc-space. 
On the contrary, the antiferromagnetic state induced by 
the effective staggered magnetic field hst along the z- 
direction can be interpreted as an ionic insulator, where 
the JW fermions stay on every other site. The conduc- 
tion band and the valence band touch at fc = only 
when (A, /ist) is taken at the origin in this A-/ist plane. 
Because of the periodicity of ^ along the fc-axis, the dou- 
ble degeneracy point at (fc. A, hst) = 0, 0) is identical 
to that of (fc, A,/ist) = (-^,0,0). 

Elementary analyses show that this double degener- 
acy point becomes the source (sink) of the vector field 
B+i (B_i) defined in the generalized crystal momentum 
space {k- A-hst space); 

BniK) = Vj^X AniK), 

An{K) = ^(n(X)|V^|n(X)), 

where K = {k,A,hst) and = {dk,dA,dh^J. Here 

\n{K)) represents the periodic part of the Bloch function 
for the valence band {n = +1) and the conduction band 
(n = —1). The vector field An,ij, changes by ^'^k^4> un- 
der the U(l) gauge transformation of these Bloch wave- 
functions; \n{K)) cxp[i(t){K)]\n{K)) , while its rota- 
tion, i.e., Bn remains invariant. Because of this gauge 
invariance, we often call the latter vector field B„ as 
a fictitious magnetic field or flux. Correspondingly, the 
double degeneracy point located at (±^,0,0) will be 
referred to as the fictitious magnetic charge (magnetic 
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monopole in the generalized momentum space), whose 
magnetic unit can be shown to be 1: 



Isi 



dS-B±i=±l. 



(7) 



Here Si represents the arbitrary closed surface which 
encloses the double degeneracy point at (fc, A, hst) = 
(+^,0,0) ^(-^,0,0). 
Based on these observations, let us consider the 

adiabatic process where the two parameters (A./ist) 
are changed along a loop Fioop enclosing the origin 
((A,/i,t) = (0,0)); 

{hst , A) = i?(cos (fi, sin (p) (8) 
R^O, ip:0^2Tr. (9) 

Then, according to the original idea of the QAPT,^^'^^'^^ 

the total number of JW fermions (/) which arc trans- 
ported from one side of this system to the other (in the 
positive direction) during this adiabatic process is equal 
to the total flux for the valence band (n = 1), B+i, which 
penetrates the 2D closed sphere spanned by Fioop and the 
Brillouin zone : 



dS'-B+i = +1. 



(10) 



loop ' 



The physical meaning of this quantized fcrmion trans- 
port is nothing but the quantized spin transport, since 
the JW fermion density is related to the ^;-component 
of the spin density. In other words, the total around 
one end of this system decreases by 1 while that of the 
other end increases by 1 during this adiabatic cycle. This 
quantization is topologically protected against the other 
perturbations as long as the gap along the loop remains 
finite, ^^'^^ in other words, as far as the double degeneracy 
points do not get out of (enter into) the 2D closed surface 
rioop X Then, we naturally expect that this 

quantized spin transport is stable against the weak ex- 
change interactions between the z-components of spins; 



Hz = |7lE[-^+(- 



.l)^A]5|5|+i. 



(11) 



In the following, we will prove that this is indeed the 
case, by introducing the exchange interactions between 
S^. As a first step, we will treat -ffxY term as a non- 
perturbed term and review which perturbations are rele- 
vant among Hz, i?dim and Hst by using the bosonization 
technique. Namely, we first introduce the slowly varying 
fields, R{x), L(x): 



fj 
R{xj) 



fee 



i{k—kF)xj 



|fe-fei?|<c«- 



,i(fe+feF)Xj 



\k+kF\<&a-^ 

Then we rewrite the spin operator in terms of these 



fermion fields R{xj) and L(xj): 



S: 



3 JO 2 

:^:R\x,)R{x,): + :L\xj)L{xj): 
+{-iy[R\x^)L{x^) + L\x^)R{x^)\, 

^t^i+i + ^3 ^j+i = fjfj+i + f]+ifj 

~ ia • [: R^{xj)dxR{xj) : — : L^{xj)dxL{xj) 



-2i{-iy 



R\xj)L{xj) - L\xj)R{xj) 



(12) 

-H.c] 

(13) 



where : ... : stands for the normal order. ''^ Here we re- 
tain the lowest order terms with respect to the lattice 
constant a both for the uniform and staggered compo- 
nents. According to the bosonization recipe, '^^ we can 
rewrite these fermion operators by using the phase opera- 
tor 0+{x) and its canonical conjugate field n(a;). Namely, 
the spin operators in eqs.(12) and (13) read; 

S^ = ^^^-{-iy^sme^{x,), (14) 



^3 "^i+i 



= —a 



^{dj+ixj))\4nn{xjf 



-{—ly — cos 9+{xj 



TTa 



(15) 



Then we substitute these equations into eqs.(4)— (6) 
and (11) and obtain the following expressions for 
HxY, Hz, Hdim and Hst in the continuum limit, where 
we neglect the rapid varying components such as 
J2ji~^y cos9+{xj) and etc.: 



H 



XY 



Hz 



h 

M [ 

a J 



2nJn{xf + ^{dj+{x)f 



(16) 



da; 



^iOMx)r + ^cos2e,ix) 



A 



(7ra)' 



■ cos 0+ {x) 



Tra- 



I 



dx cos6+{x) 



Hsi 



— ^ / dxsm0^{x). 



TTa^ 



(17) 

(18) 
(19) 



Here the origin of the third term of cq.(17) were dis- 
cussed elsewhere, '^^ which we will omit by redefining A 
in eq. (18) in the followings. Then our final expression 
for Hamiltonian reads, 



H 



dx< V 



sin(^+ + <p) + cos 20+ ] , (20) 
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where the velocity v is given hy v = {1 + ^^)- We 
also introduced the quantum parameter r] as 

which measures the strength of the quantum fluctua- 
tion. Namely, when I7I in eq.(ll) grows from 0, this pa- 
rameter decreases monotonically from 2. Even though 
eqs. (16) — (21) were derived in the weak coupling limit 
(I7I <C 1), the final form of eq. (20) is known to be valid 
until I7I reaches 1 (isotropic Heisenberg model), where 
the SU(2) rotational symmetry of correlation functions 
((5|5f) ~ (S^Si)) requires that ry = 1 at |7| = l.^e 
By using this quantum parameter r], the renormaliza- 
tion group (RG) eigenvalue of sin(6'+ + Lp) is represented 
by 2 — I while that of cos2^_|_ is given by 2 — 2r]. This 
indicates that, as long as I7I < 1 (1 < < 2) , the 
exchange coupling between the z-componcnt of spins 
{2^2^ cos2^_|_) is always irrelevant in the sense of renor- 
malization group analyses, while the dimerizing field A 
and staggered field hgt sin(^+ + ip)) are equally rel- 
evant and lock the phase operator 6^ on ^ — ip + 2iTn. 

As the system is always locked by sin(6'+ -|- (p) for 
I7I < 1, we naturally expect that the quantized spin 
transport in the case of I7I = could be generalized 
into the case of finite I7I, at least up to I7I = 1. This 
expectation is easily verified when we notice the physi- 
cal meaning of the phase operator Since the spatial 
derivative of the phase operator corresponds to the z- 
component of spin density, this phase operator is nothing 
but minus of the spatial polarization of the 2;-component 
of spins, i.e., — Ps^ = ^TfY^^=i3^j ■ The equivalence 
between these two quantities is discussed in detail in the 
appendix. Then, through the adiabatic process eqs. (8) 
and(9), (^+) decreases monotonically and acquires — 27r 
holonomy after one cycle (see Fig. 1). In other words, 
Ps« increases by 1 per one cycle, i.e., 

5Ps^ = (f dPs^ =-^f dA • dxih) = 1- (22) 

This relation always holds as far as the system is locked 
by the sliding potential sin(6'+ + <p(t)), which is true for 
111 <1. 

Then we will argue the physical consequence of eq. 
(22). Generally speaking, when the bulk system has a 
finite spin gap, the effects of boundaries range over its 
magnetic correlation length from the both ends. Then we 
naturally divide the contribution to Pg^ into following 
two parts; 

p„ _ pin , pedge pin(edge) _ 1 . , Azx 

where the ^'-summation in Pjj" arc restricted within the 
interior of the systems (r^in), while spins within the edge 
region (Oedge) contribute to P|^^**. The "interior" of the 
system is defined as the region where the spin state is 
same as that of the periodic boundary condition (p.b.c). 
Namely, the presence of the boundaries does not influ- 
ence on the spin state of the interior. On the other hand. 




Fig. 1. A sine curve represents the locking potential which slides 
adiabatically. As long as the sliding speed is low enough, the 
system (shaded circle) stays the same valley and does not jump 
into its neighboring valleys. 

the spin state within the "edge" region is strongly af- 
fected by the boundaries. Then, by deflnition, the spins 
within the interior would turn back to the same spin con- 
figuration as that of the initial state after an adiabatic 
evolution along Fioop, 

5P'^. = 0. 

Therefore the difference of the spatial polarization of 
(Sj) should be attributed to the change of spin configu- 
rations in the edge region, 

^pi'^''' = ^ E i-'5(^;) = i- (23) 

J G fledge 

In eq.(23), we next approximate j for sites around the 
right end by N/2 and those of left by -N/2; 

-\ E -5(51) + ^ E ^(^i ) - 1' (24) 

J— 2 J— 2 

where we take the origin of site index j at the center of 
the system. This approximation is allowed, since the edge 

region ranges over the magnetic correlation length from 
both ends, which is at most several sites in gapped spin 
systems. Thereby the semi-equal sign in eq.(24) could be 
safely replaced by the equal sign in the thermodynamic 
limit. When we bear in mind that the system has been 
staying on the eigenspace of S^^^^ = during this process; 

E '5(^1) + E '^(^1) = o> 

J— 2 J— 2 

eq.(24) indicates that the total around the right end 
[j ~ N/2) increases by 1 while that of the left end (j ~ 
—N/2) decreases by 1 along this cyclic evolution: 

E '^(^l) = -i' E = 

jGleft end right end 
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Namely, the spin current quantized to 1 flows from the 
left end of the system to the right end during this one 
cycle. 

In summary, the quantized spin transport is always 
assured as long as |7| < 1. This story does not alter even 
with the finite sliding speed c = as far as it is slower 
than the height of this locking potential, in other words, 
the spin gap along the loop. 

3. Numerical Analyses 

In order to confirm the above argument, we performed 
the numerical calculation in the case of I7I = 1. Namely, 
we generated numerically the temporal evolution of the 
ground state wavefunction under the following time- 
dependent Hamiltonian: 

|<^(i))=T{exp[i£di'iJ(0]}-|5)i, 

N 
i=l 

N 

+ /ist(i)E(-i)'^i' 

i=l 

27ri 27ri 
{h,t{t), Ait)) = i?(cos — , sin — ), (25) 

where T represents the time order operator, T is the 
cycle period during which the system has swept the loop 
(rioop) one time and \g)i is the ground state wavefunction 
at (/ist,A) = (i?,0). 

In Fig. 2, we show the expectation value of Sj taken 
over the final state, \(j){t = T)), both with the p.b.c. and 
with the open boundary condition. The parameter R is 
fixed to be 0.3 and J is taken as 1.5 (0.6) in Fig. 2a(b). 
The cycle period T is taken as 40 for the former and 80 for 
the latter, both of which are sufficiently long compared 
with the inverse of the spin gap observed along the loop. 

For the p.b.c, as far as the sweeping velocity is low 
enough, the final state gives completely same configu- 
rations of (Sj) as that of the initial antiferromagnetic 
ground state. However in the case of the o.b.c, the spins 
around both boundaries are clearly modified as seen in 
Fig. 2. When we read the total around the left end 
( E'=i(^|) ) and that of the right end ( E]ti20j) ) 
from Fig. 2a(b), the former increases by 0.67 (0.9) while 
the latter decreases by 0.67 (0.9) after this cyclic process. 
This result is qualitatively consistent with our preceding 
arguments. 

In order to understand the difference between the re- 
sult of the o.b.c. and that of the p.b.c, we also calculated 
the instantaneous eigcncncrgy as a function oi ip = ^ 
in Fig. 3; whore wc take the offset as the ground state 
energy. For the p.b.c. (inset of Fig. 3a) , there is always 
a finite energy gap from the ground state for all (p. On 
the contrary, the gap for the o.b.c. reduces strongly at 
(f = ^ (Figs. 3a and 3b). Furthermore, this reduced en- 
ergy gap Ai becomes smaller and smaller when we take 
the system size A'' larger (inset of Fig. 3b, N=6,8,10,12), 
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Fig. 2. (a): (5|>t=T for J = 1.5, R = 0.3 and T = 40, where 

SPs^ = {Ps-}t=T - (Ps->t=o = 0.623, (b): (5|>t=T for J = 
0.6, R = 0.3 and T = 80, where SPs^ = 0.862. 

and the four states"^^ including the ground state are 
nearly degenerate a,t (p — ^ ior N = 14. 

Let us discuss the physical meanings of this quasi- 
degeneracy. At (ft-st, A) = (0, — i?), the exchange coupling 
between S2n and S2n+i (n = 1,2,...) are strengthened 
(see eq.(25)) and dimers are formed between them, i.e., 

Si — S2 = S3 — ■ ■ ■ — Sn~2 = Sn-1 — Sn 

where Si = Si+i represents the singlet bond. Then, in 
the case of the o.b.c, the spin at j = 1 and that oi j = N 
cannot form a spin singlet due to the absence of Si ■ Sn 
term. Namely, in the limit oi J = R, the following four 
states would be degenerated, 

I i)i «) (dimer chain) ® I t)Ar, (26) 

I t)i «) (dimer chain) (g) I i)Ar, (27) 

I i)i (g) (dimer chain) (g U)Ar, (28) 

I t)i O (dimer chain) (g) I t)Ar. (29) 

The last two states ((28), (29)) do not belong to the 
eigenspace of S^^^ = and we ignore them henceforth. 
When we introduce finite J — R > 0, the ^-th order 

perturbation in terms of S1S2, S^S^ and S^-iS^ 
(and their Hermite conjugate) lifts this degeneracy. Then 
the resulting gap should be scaled as exp{—Na/l), where 
I is expected to be a magnetic correlation length. In fact, 
we can fit the size dependence of Ai by this exponential 
function as in the inset of Fig. 3b, from which the corre- 
lation length at (/5 = ^ is estimated around 3 sites. Be- 
cause of this size dependence, we expect that Ai would 
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Fig. 3. (a): Energy levels of all the eigenstates as a function of 
up = ^ for the o.b.c. (inset: p.b.e.). J = 1.5, R = 0.3 and 
N = 14. (b): Energy levels of those eigenstates with S^^^ = for 
the o.b.c. J = 1.5, R = 0.3 and N = 14. (inset: Size-dependenee 
of the minimum energy gaps Ai = miui^ [E„=i (i^) — Eg((^)] 
and A2 = min,^ [B„=2 ("p) — En—\{ip)\. Fitting curves for Ai 
and A2 are 0.32exp(-(W - 6)/3) and 0.925 exp(-(Ar - 6)/4.3) 
respectively.) 



Fig. 4. (a): Projected weight of \4>(t)) onto each instantaneous 
eigenstate |(n|(^(t))p for n = 0,1,2,4 where the w.f. |<^(t)} is 
always almost perpendicular to the excited state with n = 3: 
1(71 = 3|aHi))P < 10"^. J = l.b,R = 0.3 and TV = 14. (b): 
Traiisitiou probability P„=o^n=l as a function of cycle period 
T. Every data points are fitted by the Landau-Zencr function 

exp(— gx27r/r ) ' where c is estimated around 2.5. J = 1.5 and 
R = 0.3. 



vanish when we took the system size sufficiently large 
compared with this correlation length. 

Because of this crossing character of the energy spec- 
trum at (/3 = the system which has resided on the 
ground state (n = 0) transits, at = into the first 
excited state with n = 1 which is represented by the bold 
solid line in Fig. 3b. We can see this transition in Fig. 4a, 
whore the projected weights of \4>{t)) onto each instan- 
taneous eigenstate are given as a function of = 
Futhermore, its transition probability, P„=o^n=i, is well 
fitted by the Landau-Zener formula^^'^^ 

P_o^„=i=exp(--^^), 

as in Fig. 4b, where we plotted Pn=o^n=i for various 
cycle period T and various system size. Prom this fitting, 
c in the Landau-Zener function is estimated around 2.6 ± 
0.2.42 

As shown in the Fig. 4a, after the first transition took 

place at Lp = the second transition from the first 
excited state to the second excited state with n = 2 
(represented by the bold broken line) happens at <^ ~ 5.9. 
However its transition probability P„=i^„=2 is around 
70% and not so perfect compared with P„=o^n=i- This 
is mainly because the gap A2 around <^ ~ 5.9 between 
the first excited state and the second excited state is 
a substantial amount compared with Ai (see Fig. 3b). 



However, as in the inset of Fig. 3b, its size dependence 
is also scaled by 

A2 - exp(-7Va/0, (30) 

where the magnetic correlation length I is estimated 
around 4.6 sites. Therefore, when we take the system 
size large enough, wc can naturally expect that the gap 
A2 would reduce exponentially, which enhances the tran- 
sition probability P„=i^„=2 drastically. 

In order to verify this expectation, we perform another 
numerical calculations. As it is very difficult to calcu- 
late numerically with larger system size {N > 14), we, 
instead, decrease the ratio J/R in order to reduce the 
magnetic correlation length I, which is also expected to 
reduce these gaps according to eq.(30). The result for 
J = 0.6, R = 0.3 and = 14 is summarized in Fig. 
5. FigureSa indicates the energy spectrum for the lower 
four excited states with S^^^ = and Fig. 5b shows the 
projected weight, |(n|(/)(t))|2 for n = 0,1,2,4. Figure5c 
represents the expectation value of Sj taken over the 
|n = 4)(p=27r=o as a function of site index j. Then, as is 
expected, the gaps we observed in Fig. 3 (Ai,A2) de- 
crease drastically. In addition to them, there appears a 
new gap A3 between the excited state with n = 2 and 
that with n = 4 a.t Lp = 6.15. This gap A3 is also very 
tiny and scales as exp(— A^/2.25). 

As a result of these crossing characters of energy spec- 
trum a,t ip = ^,5.97 and 6.15, three Landau-Zener type 
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Fig. 5. J = 0.6, R = 0.3 and N = 14. (a): Energy spec- 
trum for the cigenstates with S^^^ = as a function of ¥',(b): 
Projected weiglit of |</>(i)} onto the lower five eigenstates, i.e., 
|(</i(t)|n)p(n = 0, 1,2,4), (c): Expectation value of 5| taJsen over 
the fourth excited state with n = 4 |n = 4)ip=27r 



tunnelings take place as shown in Fig. 5b. Namely, the 

wavcfunction \4>{t)) changes its weight from the ground 
state with n = to the first excited state with n = 1 
(represented by the bold solid line in Fig. 5a) at <^ = ^ 
, from the first excited state to the excited state with 
n = 2 (bold broken line) at = 5.97 and lastly from the 
excited state with n = 2 to the fourth excited state with 
n = 4 (bold dotted line) at Lp = 6.15. All these tran- 
sition probabilities are almost 100% in accordance with 
our previous expectations.^^ 

After all, through these transitions, the wavcfunction 
is raised from the ground state onto the fourth excited 
state (n = 4) during this last quarter of the cycle. When 
wc see the expectation value of Sj taken over this fourth 
excited state at <^ = = 27r (see Fig. 5c), it is almost 
identical to the solid line of the lower panel in Fig. 2, 
where the total on the left end {J2^^i{Sj)) decreases 

by 0.9 while that of the right end {J2]=i2i'^j)) increases 
by 0.9 in comparison with that of the initial ground state 



wavcfunction. 

To summarize, the numerical observations found in 
Fig. 2 can be ascribed to the difference between the en- 
ergy spectrum with the p.b.c. and that with the o.b.c. In 
the latter certain excited state decreases its energy 

level until ip = ^ and then picks up the system which 
has been staying on the ground state. Then, through 
the Landau-Zener tunnelings, several excited states carry 
the system hand in hand perfectly upon a particular ex- 
cited state at = 0, whose Ps' increases by almost 1 
in comparison with that of the initial state. These ob- 
servations remain unchanged if the sweeping velocity is 



always higher than 

— > max . 

T i=l,2,3V c / 



(31) 



This lower limit for the sweeping velocity is expected to 
vanish in the thermodynamic limit, since the Ai^2,3 re- 
duces exponentially to zero as the system size A'' becomes 
larger. Therefore the conclusion in cq. (22) is consistent 
with the numerical results of this section. 

Finally let us mention about the relevance of our story 
to real systems such as Cu-benzoate and Yb4As3. There 
the effective staggered field along the z-direction ac- 
companies the uniform effective Zeeman fields along the 
a;(or y)-dircction, HZ^^\ Furthermore, since these sys- 
tems break the bond-centered inversion symmetry, there 
is a Dzyaloshinskii-Moriya (DM) interaction. Its DM vec- 
tor must alternate bond by bond in the absence of the 
electric field. In addition to them, the uniform compo- 
nent of the DM vector would be also induced by the 
applied electric field in general. Even if these terms arc 
gradually introduced to our model calculations, 69+ re- 
mains invariant as far as the gap along the loop remains 
open, in other words, the system along the loop remains 
locked by the sliding potential sin(^+ -|- ip{t)). This is a 
manifestation of the fact that the quantized 56+ (in the 
thermodynamic limit) is a topologically protected quan- 
tity. However, detail studies on the effect of these terms 
might be remaining interesting problems. 

When the system travels around this critical ground 
state M times, total z-component of spin around one end 
increases by M, while that on the other end decreases 
by a same amount. Such an inhomogcneous magnetic 
structure can be detectable around sample boundaries 
and/or magnetic domain boundaries by using some op- 
tical probes. However, we must mention that the excited 
state with Ps^ = m (In the case of m = 1, this corre- 
sponds to |n = 4)^=0 in Fig- 5) might fall into the ground 
state \n = O)ip=o by way of the .x-componcnt of the uni- 
form magnetic field at around t — niT, when we 
sweep the system slower than 2hbH^- In order to avoid 
this relaxation process, the sweeping velocity should be 
taken faster than 2ijlbH^- 



— > 2hbH^, 



(32) 



However, in order to keep the system on a particular min- 
imum of the locking potential cos(6'+ -|- ^), the velocity 
must be also slower than the spin gap observed along the 
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loop, 



(33) 



Therefore we need to have a finite window between 
the lower limit of the sweeping velocity and the up- 
per limit. In the case of Cu-bcnzoatc, the spin gap of 
O.lmeV is induced by the applied magnetic field of IT 
(2/Us X IT ~ O.lmeV). Thereby, the lower limit of this 
sweeping velocity is imhappily comparable to the higher 
limit in this system. However, spin gaps around O.lmeV 
could be also induced by the magnetic field smaller than 
IT in those quantum spin chains with relatively larger 
J, where substantial amount of the window between the 
upper limit and the lower limit would exist. Therefore, 
it is not so hard to realize our story experimentally on 
the quasi-lD quantum spin systems with relatively large 
intrachain interaction and with two crystallographically 
inequivalent sites in its unit cell. 

We also want to mention about the magnitude of the 
electric field required in order to induce the dimerization 
gap enough to be observable and also enough to be com- 
parable to the spin gap induced by the magnetic field, 
which is around O.lmeV. Since it goes beyond the scope 
of this paper to quantify microscopically the change of 
the exchange interaction induced by external electric 
fields, we will pick up some reference data from other 
materials. Let's see magneto-electric materials, where the 
applied electric field often changes its effective exchange 
interactions J due to its peculiar crystal structure and its 
change A could be quantitatively estimated via the re- 
sulting magnetization.*^ In the case of famous magneto- 
electric material Cr203, the intra-sublatticc exchange in- 
teraction for the + sublattice and that for the — sublat- 
tice should be equal to each other because of the inversion 
symmetry ( J+ = J_ = J) . Then an electric field applied 
along its principle axis breaks this symmetry and induce 
a finite difference A = J+ — J_, where A/ J was esti- 
mated around 10~^ under E ~ IkV/cm.*^ Namely, 

J = cE. (34) 

with its linear coefficient 10~^[cm/kV]. When we apply 
this coefficient to our Heisenberg model, eq.(2) with J ~ 
O.leV, we need an electric field of order of IkV/cm in 
order to induce the dimerizing gap J^/'^A^/'^ of order of 
O.lmeV, which is still one order of magnitude smaller 
than the typical value of the Zener's breakdown field (> 
lO^V/cm at T < lOOK) of ID Mott insulators.'**' 

4. Conclusion and discussion 

In this paper, we propose the quantum spin pump- 
ing in S = 1/2 quantum spin chains with two crys- 
tallographically inequivalent sublatticcs in its unit cell 
(T, /bond '■ broken and /site : unbroken) and with rela- 
tively large intrachain exchange interactions. Due to its 
peculiar crystal structure, an applied electric field (E) 
and magnetic field (//) endow its critical ground state 
with the finite spin gap via the dimerizing field (A ~ E) 
and the staggered magnetic field (ft-si ^ H), respectively, 
where the phase operator 9-^- in sine-Gordon theory is 



locked on a particular valley of the relevant potential 
^/^^+^sm0+ + If) : {§+) = ^-ip + 27rn. When the 
system is deformed slowly along the loop which encloses 
the critical ground state in the E-H plane, the lock- 
ing potential + A^ sin (6'+ + ip) in the sine-Gordon 
model slides gradually (see Fig. 1). After one cycle along 
this loop, the expectation value of the phase operator 9^ 
for the ground state acquires 27r holonomy. This means 
that a spatial polarization of z-component of spins Ps^ 
increases by 1 after this cycle. In other words, the z- 
component of the spin current quantized to be 1 flows 
through the bulk system during this cycle. 

These arguments are supported by the numerical anal- 
yses, where we performed the exact diagonalization of the 
finite size system and developed the ground state wave- 
function along this loop in a thermally insulating way. 
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Appendix: Equivalence between the phase oper- 
ator and the spatial polarization op- 
erator 

In this appendix, we will show that the contin- 
uum limit of Pgz corresponds to the phase operator 

^Wa I ^^2n ^ where related work was done by Naka- 
mura and Voit.^^ Historically the polarization operator 
is known to be an ill-defined operator on the Hilbert 
space with the periodic boimdary condition (p.b.c.) and 
is consequently formidable to treat in solid state physics. 
However, recently King-Smith and Vanderbilt^^'^^ and 
Resta^*"'"'"' showed that the derivative of the polariza- 
tion with respect to some external parameter A can be 
expressed by the current operator which is in turn well- 
defined in the Hilbert space with the p.b.c. Based on 
this observation, they quantitatively estimated electronic 
contributions to macroscopic polarizations in dielectrics 
and semiconductors such as GaAs, KNbOa and III-V ni- 
tride. Following their strategies, we will show the equiv- 
alence between the derivative of —^-^ J dy(^^^M^ with 
respect to A and that of (Ps^)- 

According to the standard perturbation theory,**^ the 
derivative of the polarization with respect to some me- 
chanical parameter A is given by 



dX 



■E 



(ff(A)|[Ps.,g]|jV(A))(Ar(A)|g|P|g(A)) 
{E,{X)-E^{\)Y 



c(A X) 



where we assume the system has a finite spin gap from 
the ground state. The time derivative of Ps^ in the nu- 
merator reduces to the summation of the current opera- 
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tor over all bonds 



= 4 [i^l' [-^ + (-1)'^] i^^S^+i + ^•^•) 

j=i k=i 

. N-l 



that this is indeed the case, wo will carefully reexamine 
the bosonization of Sj'Sj]^-^ + ^^^^^"^^ : 



N N-l 



5+57^, +H.C. 



H.c. 



Thereby the contribution of every single bond is always 
0{1/N), which is not the case with Pg^ = jj- ^2^=1 J^j- 
As a result, we are allowed to cistimate with the periodic 
boundary condition the r.h.s. of eq. (A-2) in the thermo- 
dynamic limit, which we should have figured out with 
the open boundary. This is because their difference is 
attributed to the spin currents on the several bonds^'' 
around the boundaries, which is at most 0{1/N) in our 
spin gapped system. Then we will consider eq. (A-2) with 
the p.b.c. and rewrite it in terms of the bosonization lan- 
guage. That is to say, we express the current operator in 
eq. (A-2) in terms of the phase operator 6+{x) and its 
conjugate field Ii{x), 



(r.h.s of eq.{13))-i{-iy—-[R^Xj)djL{xj) 



-Ll{xj)da:^R{xj)-'H..c. 
= (r.h.s. of eq. (15)) 



+ TT 



n {xj ) - cos 6+ {xj ) 
+ n(xj) • cos^+(a;j) ■ il{xj) + H.c. 



(A-7) 



[R^ixj+^)e-"^^ 
X \^R{xj)e''""'^ + L{xj)e- 



-H.c. 



Then, by taking its staggered components, we obtain the 
following expression for Hdim, instead of eq.(18), 

A / dx\ ^-:tCos^+ — (dxO+] cos^+ 

J 7ra^ Stt V / 



-2i- 



R^{xj)R{xj)-L^{xj)L{xj) 
{-iy\R^{xj)dxL{xj) - L\xj)dxR{xj) + H.c 

= 4m{xj) 

+ i(-l)^[2fl(xj)cos^+(xj) +H.C. 



- cos 61+ + n - cos 6*+ • n + h.c. 



.(A- 



Accordingly, the time derivative of Ps^ in the continuum 
limit can be expressed in the following form. 



Here wc want to mention that, irrespective of this modi- 
fication, cos 6+ in iJdim always locks the phase operator 
in combination with sin^+ in hgt as far as |7| < 1 (see 
the text). This is because the RG eigenvalues of the addi- 
(A-3) tional terms such as (dxO+Y cos^+, (H)^ cos^+ and etc. 
are all negative and thereby irrelevant in the sense of the 
renormalization group study. However these additional 
terms in ifdim cannot be discarded, since the commuta- 
tor between these terms and the phase operator produces 
the last two term of eq.(A-4): 



2jn-|- A 



ncos^+ -FH.c. 



k.4) 



,^0+{x) 



Prom now on, we will show that this expression is iden- 
tical to the time derivative of the phase operator: 



2ti 



= -A 



n(a;) cos 6+ {x) + cos 6+ {x)Il{x) 



(A-9) 



(r.h.s. of eq.(A-4)) = i[^ 



dx^,H]. 



(A-5) 



The first term in the r.h.s. of cq.(A-4) comes from the 
commutator between 6+ and Hxy given in eq.(16). 



2n 



XY 



After all, by comparing eqs.(A-6) and (A-9) 

with cq.(A-4), wc can safely replace —i[Ps^,H] by 

■^Jdx[^,H] in the continuum limit and rewrite 
eq. (A-1) into a following simple form by using the 
bosonization language: 



= -2jn(a;). 



(A-6) 



On the other hand, the commutator between the phase 
operator and i?dim given in eq.(18) vanishes and does 
not produce the last two terms of eq. (A-4). This is be- 
cause cq.(18) is not an accurate expression for -ffdim in 
the continuum limit and needs some additional terms, 
whose commutators with the phase operator correctly 
yield the last two terms of eq. (A-4). In order to prove 



dx 



Na 



Njtg 



{g\['-^,H]\N){N\^\g) 



{Eg - 



+ c.c. 



{g\^-^\N){N\^\g} + C.C. 



d , 0+ix) , 

27r ' 



(A-10) 
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